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A symmetric binary mixture (A,B) below its critical temperature T c of unmixing is 
considered in a thin film geometry confined between two parallel walls, where it is as- 
sumed that one wall prefers A and the other wall prefers B. Then an interface between 
the coexisting unmixed phases is stabilized, which (above the wetting transition tem- 
perature) occurs in the center of the film for an average concentration of c = 1/2. The 
problem is considered how the concentration profile c(z) across the thin film depends on 
the film thickness D. By Monte Carlo simulation of a lattice model for a polymer mixture 
it is shown that for relatively small D the width of the interface scales like w oc D while 
for larger D a crossover to a behavior w oc \/Z) occurs. This behavior is explained by 
phenomenological theories: it is shown that the behavior at small D can be understood 
by a suitable extension of the Cahn Hilliard "gradient-square" -type theory, while the 
behavior for large D can be traced back to the behavior of capillary waves exposed to 
a short range potential by the walls. Corrections due to fast concentration variations, 
as they occur in the strong segregation limit of a polymer mixture, can be accounted 
for by self-consistent field theory. Subtle problems occur, however, with respect to the 
proper combination of these theories with the capillary wave approximation, particularly 
at intermediate values of D. 



I. INTRODUCTION 



The theoretical understanding of interfacial profiles between coexisting phases has been a longstanding 
challenge JT|-f26| . After Cahn and Hilliard Q in a seminal pioneering work proposed the " gradient square" 
theory for binary mixtures, to derive the (intrinsic |l[) interfacial profile in the framework of a mean field 
type treatment, this approach has been extended to treat interfaces in polymer solutions 0, polymer 
blends [p -[l0 Jrc| - [T^ , ^6| - p3| , p5| , p6| block copolymer mesophases and other modulated phases in complex 
fluids |4|Jl5|2J etc. Of course, for a full understanding of the interfacial profile one has to amend the 
mean field treatment by a consideration of statistical fluctuations neglected in this treatment, notably 
long wavelength fluctuations of the local position of the center of the interfacial profile (usually termed 
"capillary waves" § f|l|]2|[|^§3 ) . 

Here we focus on another extension of the " gradient square" theory, namely the profile of an interface 
in a geometry confined between two parallel walls a distance D apart. While interfaces bound to external 
walls have been considered for a long time in the context of wetting transitions P,|3|, [^ , j32| -|38fl , interfaces 
confined by two walls have been studied only more recently |39|-|46|] , focusing on the interface localization 
- derealization transition that may occur in this geometry. In the present study, however, we investigate 
another aspect of these confined interfaces that has been discovered in recent simulations |47|-|49|] and 
experiments on polymer mixtures |^7],^0| : there occurs a significant reduction of the intrinsic width wq 
of the interfacial profile already for relatively thick films, D ^> wq. We expect that this phenomenon is 
important also for a description of interface-location - derealization transitions and wetting phenomena, 
since it implies a renormalization of parameters entering the effective interface Hamiltonian | ]5l| . This 
"squeezing" of the intrinsic interfacial profile also is important for experimental studies that try to 
deduce effective interaction parameters from measured concentration profiles 52 of polymer mixtures, 
but apply the standard theory for interfacial widths of strongly segreg ated systems that 
ignores the finite thickness D of the thin film that is used. 

The outline of the present paper, hence, is as follows: we first (Sec. 2) present a phenomenological 
theory of this squeezing of interfacial profiles in thin films, by a fairly straightforward extension of the 
theory of Parry and Evans Q], who did not pay attention to this effect. As these authors, we consider 
only short range forces of the walls for simplicity, although for a quantitative description of experiments 
0,|O|J52|5| long range van der Waals forces should be used (numerical evidence shows that for long 
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range forces there even is a stronger squeezing of the interfacial profile due to confinement than for 
short range forces 0)) . We show that there is a rather broad range of D where u>o varies nearly 
linear with D. While this treatment implicitly implies the case of a weakly segregated mixture, where 
concentration variations are slow, we then (Sec. 3) treat the alternative case of a strongly segregated 
polymer mixture where the profile varies across the interface from <pA = to 4>A = 1> but nevertheless 
a "square gradient" theory applies, in the limit of long chain length of the polymers The case 

of intermediate segregation is more complicated, since even a free, unconfined interface of a polymer 
mixture exhibits a profile involving two different lengths |12j| . We study such intermediate cases in Sec. 
4 by a numerical version of the self-consistent field theory [||J|Jll|,|l|,|^-[l|| that uses Monte Carlo - 
generated single chain configurations as an input [|T]j5^j5J]. Sec. 5 then summarizes and discusses our 
results, paying attention also to the problem how the present results can be combined with the effects 
of capillary waves in this confined geometry (which for D — * oo lead to an apparent j|7 48 1 interfacial 
width w oc \f~D while the intrinsic width wq saturates at a finite value - the mean field result for an 
unconfined free interface). 



II. REDUCTION OF INTRINSIC INTERFACIAL WIDTHS OF COEXISTING 
SYMMETRICAL MIXTURES CONFINED BY "COMPETING" WALLS 

We consider here a symmetrical binary mixture (AB) sufficiently close to the bulk critical point such 
that the concentration variations of interest are of rather long range, and a Ginzburg-Landau type 
description in terms of the order parameter m(z) = {[c(z) — c cr it] /c cr ;t}(here c(z) is the concentration 
of species A at position z and the critical concentration c cr i t of species A is c cr i t = 1/2) is applicable. 
We consider a thin film geometry with two "competing" planar walls, oriented perpendicular to the z- 
direction, a distance D apart (Fig. 1). "Competing" walls mean that the left wall preferentially attracts 
species B, and the right wall preferentially attracts A, and thus the concentration c(z), as well as m(z), 
increases monotonically from left to right. We take these forces between the walls and the molecules (or 
atoms, respectively) of the mixture of short range, i. e. they act locally on particles adjacent to the walls 
only, for the sake of simplicity, and furthermore we assume their absolute strength equal, so that the 
profiles m(z) shown in Fig. 1 must be antisymmetric around the point z=0, m(z=0)=0 i.e. the center 
of the profile in the middle of the thin film. Of course, here we have also assumed that these forces 
exerted on the particles on the wall are sufficiently strong, so that no symmetry breaking due to the 
localization of the interface is possible at the considered temperature. Note that for large enough D the 
interface localization-delocalization transition would occur close to the wetting transition temperature 



of the corresponding semi-infinite system J39-46|, and recall the argument due to Cahn |32| that close 
to the bulk critical point of a mixture the surfaces always should be wet. 

We now emphasize that in general we expect five distinct regions in m(z) given the assumptions made 
above: In the first region close to the left wall, m(z) will decay from its value —mo(D) right at the 

wall to the value — m&, characteristic for the left branch of the coexistence curve, Ccoex, of the binary 
mixture in the bulk |—TOfc = Ccocx — c cr it /c cr it, +mf,= Ccocx — c cr it /c cr it|. This decay occurs over 
a length scale of the same order as the correlation length in the bulk. In the second regime where 
m(z) ps — rrifc {centered at —Zb defined by m(—Zb) = —mi,} the slope of the profile dm(z)/dz reaches 
a minimum, and in the limit D— > oo the profile would become completely flat. Near the center of 
the film we have the interfacial profile in a strict sense, which for large D is described simply by Jl|-^| 
m(z) — nib tanh [z/wq(D)] with wq(D — > oo) = 2£&. For not so large D, however, the "intrinsic width" 
wq(D) of this profile is significantly reduced, cf. Fig. 1 (lower part), and this interfacial regime is no 
longer so clearly distinguished from the first regime, since the second regime where m(z) stays nearly 
flat when m(z)« — nib, nearly has disappeared {and so does the fourth regime where m(z)« +nib near 
z=Zb}. In the fifth regime we have the increase from m(z)« +rrib to m{z = D/2) = mo(-D), which clearly 
is dependent on the precise nature of the boundary conditions at the walls (unlike nib and which are 
not dependent on these boundary conditions, of course. 

From this (qualitative) expectation that the concentration profile exhibits five distinct regimes (the 
central regime is expected to have a width of 2£f,, the others must have at least a width of (j> in order for 
this picture to make sense) one already can predict that the asymptotic regime of this interfacial behavior 
is only reached for D 3> 6£f,, while a rather nontrivial behavior {including a significant reduction of the 
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width Wo(D) in comparison to 2^} is expected if this condition is not met. Of course, Fig. 1 makes 
sense only for the case of rather weak segregation, close enough to the bulk critical point, where m(,=l 
for our normalization). If one can reach a case of strong segregation but still is in a regime of complete 
wetting, then nib — * 1 and also too (D) — 1, and one obtains a much simpler situation than shown in 
Fig. 1, with m(z) varying from —1 at z=—D/2 to +1 at z = D/2 in a single step {a profile varying 
similar to tanh [z/wq(D)\\. It has been predicted |35| that such a situation (strong segregation in the 
bulk but complete wetting of a surface) is typically realized in polymer mixtures, and this prediction 
was confirmed by corresponding Monte Carlo simulations jl8],[36| . We shall return to this case of strong 
segregation in the next section. 

In order to calculate the explicit form of this nontrivial profile in Fig. 1 from a Ginzburg-Landau type 
theory p]-|6|j42]] , it is convenient to rescale the order parameter m(z) by its bulk value, M(Z) = m(z)/mb, 
and also lengths are measured in units of the bulk correlation length, Z = z/£b. Then the rescaled free 
energy of the film can be written as follows [p|,ff2"[ 

-^M 2 (Z) + ^M\Z) + {dM/dZf 

-D/{2( b ) 

+ (6/A) [M 2 (-D/2) + M 2 (D/2)] + Hi [M(-D/2) - M(D/2)] (1) 

Note that at z — +D/2 we have applied a (normalized) field Hi and at z = —D/2 a (normalized) field 
—H± coupling to the order parameter; this term represents the preferential attraction of species A to the 
right wall, species B to the left wall. As is well known |§, one must allow in the "bare" wall free energy 
{the second line of Eq. ([!])} a term proportional to iVP as well (to account for the effect of "missing 
neighbors" and possible change of interactions near the wall, etc.). Following common notation 0|, the 
coefficient of this term has been put inversely proportional to the so called "extrapolation length" A which 
is a microscopic length (of the order of the interaction range in a model with short range interactions). 
Note that F(D) is normalized per unit surface area of the walls, and in units where the bulk free energy 
per unit volume is —1/4. 

Now the profile M(Z) in Fig. 1 results from seeking the minimum of the free energy functional, Eq.(jl]). 
The Euler-Lagrange equation corresponding to Eq.ffl) is 

- M(Z) + M 3 (Z) - 2 d2 fl ( ? Z) = (2) 
dZ z 

with the boundary conditions 

dM/dZ- (&/A)M = Hi/2, Z = -D/2 , (3) 



^F(D) = J dZ 



dM/dZ + (&/A)M = Hi/2, Z = +D/2 . (4) 
After multiplication of the Euler-Lagrange equation with M = dM/dZ one can integrate once to find 

{if) 2 = -\ M2{Z) + \ M4(Z) + \- A P( D )= [M 2 ( Z ) - 1] 2 /4 - Ap (D) , (5) 

where the integration constant was denoted as 1/4 — Ap (D) , for the sake of consistency of notation 
with the work of Parry and Evans p2| . For D — > oo, we expect that dM/dZ — > for M(Z) = ±1 
(corresponding to the flat regions 2 and 4 in Fig. 1) and thus Ap(D — > oo) — > 0. 

Eq.(||) can also be combined with the boundary conditions, Eqs.(||), (^), and hence we conclude, using 
M = M(Z = D/2&) as an abbreviation, 

(M&/A - ffi/2) 2 = [M - l) 2 /4 - Ap(D), (6) 

and from Eq.(||) we find, invoking the symmetry M(—Z) = —M(Z), 

M 

D/2£ b = J dM/y/(M*- l) 2 /4 -Ap. (7) 
o 
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Using once more Eq. (g) one finds a closed expression for M as function of D in terms of the inverse 
function expressed as a quadrature, 



M 



= / dM/J (Af 2 - If - (Af 2 - l) 2 /4 + (M&/A - H x /2) 



(8) 



The maximum slope of M(Z) at Z=0 {where also M(Z=0)=0, see Fig. 1} can be written as 



dM 
~~dZ 



= [1/4 - Ap{D)] 1/2 = (M&/A - i7i/2) 2 - M 4 /4 + M 2 /2 



1/2 



(9) 



Noting that for D — > oo we must have Ap = (otherwise the integral in Eq. (j7|) would converge to 
a finite constant) we can obtain from Eq.^) straightforwardly the local order parameter at the 
surface, for this situation of complete wetting, 



M006/A - Hi/2 = ~(Ml - l)/2. 



(10) 



The (second-order) wetting transition occurs when j55| Moo — > 1, i. e. for Hi c = 2£ b /\. Hence the 
situation shown in Fig. 1 requires H± > H\ c and then M^ > 1. In general, we expect for finite D that 
Af < Mao and for 00 we have a smooth convergence of M towards Moo- This expectation can be 
verified from Eq.(o) by writing M = — SM and expanding to first order in SM, using also Eq.(10), 



2 (£&/A) Moofo/A -Hi/2)+ M oc (M4 - 1) SM = Ap, 



(11) 



SM = 



{-S, b /X + M x ){Ml 



1)1 'Ap. 



(12) 



For large D we shall find that Ap is very small {Ap ot exp(— D/2£b), see Eq.([T8|) below}, and thus 
Af as Moo. For we see from Eq.(@) that then also M — ► 0, and hence in this limit the integral can 
be evaluated expanding the square root as 



Af 4 M 2 1 



-1/2 



Z -Ap 



-1/2 



[l + Af 2 /(1 - 4Ap) -+....] 



and hence we obtain 



M/{--Ap 



1/2' 



In this limit where D/2£ b w Af / (± - Ap) 1/2 < 1 we also have M(Z) 
simply linear, and from Eq.(jJ) we conclude that 



1 + Af 2 /(3- 12Ap) 



(13) 



(14) 



d M , /, AN 

— | m ax= (l/4-Ap) 



1/2 



Hi/2, 



D -> 0. 



2MZ&/D, the profile is 



(15) 



We expect that Eq.(|l5[) remains true as long as Af <C 1, i. e. D/2£ b <C 2/ 



Ap«4/ffi. Note that 



in our normalization Hi is dimensionless and typically larger than one, since ffi > ffi c = 2£&/A ^> 1. 
The fact that the inverse of dM /dZ\ max , which can be considered as a measure of the (normalized) 
interfacial width wq(D), remains non- vanishing as D — > can be understood from the fact that the 
gradient energy (dM/dZ) 2 in Eq. ([!]) disallows too steep gradients and hence for D — > one can no 
longer have a profile from —m b to +m b in Fig. 1 but only from —in to +rh where also rh = Mm b — > 
as D— > 0, because otherwise (dM/dZ) 2 would diverge in this limit. 

We now consider the inverse limit of Eq. (^) , namely D — > 00 when Af — ► . This means that the 
profile in Fig. 1 develops a very broad plateau near 771(2: ) « ttk, for z > {and m(z) ps — for z<0, 
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respectively}. As a consequence, the dominating part of the integral in Eq.(^) comes from M«l, and 
this suggests to approximate the integral as follows 

n r 1 , r AI , 

— = / dM/y/(M 2 - l) 2 /4 -Ap + / dM/^(M 2 - l) 2 /4 - Ap 
2 4& Jo Ji 



dM/y/(M 2 - l) 2 /4- Ap+ / dM/V(M - l) 2 - Ap 
o Jl 



M 



dM/ y/(M 2 - 1)2/4- Ap + In 



M-l + y(M-l) 2 -Ap 

7^ 



(16) 



The first integral in the last line of Eq.(16) is independent of the boundary condition, and tends to 
In [C/\/— Ap] for Ap — ■> 0, where a numerical evaluation shows that the constant C is roughly C«4. 



The second term in the last line of Eq.(|l6|) becomes in this limit In 2(M — l)/y/—Ap 
conclude 



26 



In 



2C(M-1) 
-Ap 



and hence we 



(17) 



-Ap= 2C(M - 1) exp(-£>/2&). 
Hence Eq.(^) yields in this limit for the maximum slope of the profile at Z = 
dM 



dZ 



1/4 + 2C(M - 1) exp(-D/26) = 1 + 8C ( M ~ X ) exp(-_D/2&) 



From Eq.(|19|) we obtain our central result for the width wq (D) in Fig. 1, namely 
w (D) 



w (oo) 



= V 1 + 8C(M - 1) exp(-£>/2& 



(18) 



(19) 



(20) 



Eqs.([19|) and ( p0| ) describe the approach to the limit of an "intrinsic" interface in the bulk, which has the 
width wo = wo (oo) = 2£t in our units. Since M can appreciably exceed one, and the constant 8C ~ 32 
is rather large, the approach of wo (D) to wo (oo) is rather slow, and hence a significant reduction of the 
interfacial width of a confined interface is predicted. Since in the considered limit we may replace M by 
Mqo in Eq.(EOj), we can use Eq.([10|) to write, for i/inear Hi c , using also C«4, 



w (D) 



w (oo) 

In Fig. 2 the approximation, Eq.(|2< 



= y/l + 16(ffi - H lc ) exp(- J D/26)/(l + H lc /2). 



(21) 

is compared to a full evaluation of Eqs. (||) and (p|) which can 
only be done numerically: for a given choice of £f,/A an d Hi we first obtain M^ from Eq.niO). Then we 
take for this choice of £&M an d Hi several values of M < M w and evaluate the values of D/2£t that 
then correspond to each M. Since dM/dZ\ max is readily given in terms of M as well, from Eq.(^), one 
can easily plot the ratio wo(D)/wq(co) = dM / dZ\ maXt £> =00 / dM / dZ\ maXi D versus D/2^. From Fig. 2 we 
see that for large -D/6 the reduction of wo (D) /wo (oo) is only dependent on M^, as expected from 
Eq.(po[) since in this limit M M^ can be used there. In addition, for 2 < D/(j> 8 the variation of 
wo (D) with D is approximately linear, and a saturation of wo (D) at wo (oo) only occurs for £>/£*> — 12 
, as predicted on the basis of the qualitative discussion of Fig. 1 at the beginning of this section. 

If for a given M the corresponding D has been evaluated from Eq.^), Ap(D) being known from 
Eq.(^|), one can use Eq.(||) to compute the full profile M (Z) in analogy with Eq.(^), namely 

-M{Z) 

dM'/y/(M' 2 - l) 2 /4 -Ap (22) 

Putting M (Z) = 1 here one obtains the value Zf, = Zb/£b defined in Fig.l. All the details of Fig. 1 
thus can be verified explicitly. 
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III. SELF-CONSISTENT FIELD TREATMENT FOR CONFINED INTERFACES OF 
STRONGLY SEGREGATED POLYMER MIXTURES 



The Monte Carlo simulations of a lattice model of a strongly segregated polymer mixture between 
competing walls J48| already have provided numerical evidence (Fig. 3) that also in this limit there 
is a regime where the observed interfacial width w apparently scales as w oc D, despite the fact that 
interfacial fluctuations yield a strong increase of w with D for D 3> wq (oo) as well. 

In order to present a theoretical understanding for this problem of interface squeezing by confining 
walls in strongly segregated mixtures, we present here an extension of the theory due to Helfand et al. 
|^|,^| The basic quantities are the probability densities qa (z, t), q B (z, t) that one end of an A or B chain 
with degree of polymerization t is at z. The position of the other chain end is arbitrary. For a blend 
of A and B homopolymers of the same chain length Na = N B = N, these functions satisfy modified 
diffusion equations, 

d b 2 d 2 q A (z,t) 

—q A (.z,t) = — — — -2 w A qA{z,t), (23) 



d , . b 2 d 2 q B (z, t) , . . ,. 

—q B (z,t) = - > -w B q B (z,t), 24 

ot 6 az A 

where the "fields" wa,wb derive as w a = df/dp a from a free energy density that depends on the 
densities p A (z), Pb{z) 



,-N 

N 



= dtq a {z,N-t)q a (z,t), a = A,B, (25) 



as 



f = XPAPB + ^(pA + PB-l) 2 - (26) 

Here the factor (k B T) 1 is absorbed in the free energy density, \ is the Flory-Huggins parameter that 
causes the unmixing of the polymer mixture (x-/V 3> 1 defines the strong segregation limit), and £ 
controls the inverse compressibility of the blend. 

While Eqs.(|2^) - (^6|) define the general framework of the self consistent field theory of polymer 
blends, we here are only interested in the limit N — > oo and treat also the blend as incompressible. As a 
consequence of the limit N — > oo, one can take |p|,[9| p a = q 2 (a = A, B) and put dq a /dt = 0. Redefining 
the units of length such that b/y/6x = 1, Eqs.([23|) and ( p4| ) then can be replaced by 

d 2 q a {z) 19/ . . 

-^--2dqrr^ a=AB - (27) 

We can interpret Eqs.(|2^) again as Euler-Lagrange equations of the Lagrangian 

r - {!«„,„, + + !(£)'}. ,2s, 

where we have a " conservation law" for the Hamiltonian 

-h(qA,q B ), (29) 



In the inhomogeneous case the quantity of interest is the excess free energy ( again units of (k B T) 1 
being used ) 

Fcxc = / dz(f - w A p A - w B p B ) (30) 



H= l - 




l( 


'dq B \ 


2 






v dz / 


as a "time" 


derivative. 


as 


usual. 
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which after integrating by parts can be rewritten as 



F exc = -2H( Zl -z ) + 2 



dg A 
dz 



dz 



qA 



dqA 
dz 



qB- 



dqi 



(31) 



We now simplify the problem by requiring strictly local incompressibility pa + Pb — q\ + Qb = 
everywhere in the system. It is convenient to express this condition using polar coordinates 



qA = sin <p, qs = cos < 
which amounts to replace the Lagrangian in Eq.(|23) by 



dz 



2 V dz 



where 



f{4>) = PaPb = - sin 2 (20). 
Then the excess free energy formulated in Eqs.([30|) and ( pl| ) becomes 

F cxc = -(zi - z )A/4 + 2^ 1 dz 

where A is a constant of motion resulting as 

A/8 



1 

2 V dz 



(32) 



(33) 



(34) 



(35) 



(36) 



This constant A is the analog of the constant —Ap(D) of the previous section. 

As a first step, we apply this formalism to an interface in an infinitely thick system without any 
boundary effects, such that \im z ->± 00 dq a /dz = 0, and also lim z ^± DO f{qA,QB) — 0. Since thus the 
constant of motion {Eq.(]2^)| H — and also A = , we have 



+ °°dz^V 



Vdz J 

and the corresponding Euler-Lagrange equation 

^ 1 • OA 

— = - sin 2m 

dz 2 Y 



7x12 dz 



(37) 



(38) 



is solved by 

ln(tan^) = z, <f> — arctan(exp(z)) (39) 
and noting Eqs.(^) this is recognized as the familiar | [l + tanh profile for the density, || 

exp(z) 



qA/qB = exp(z), 



QA 



1 + exp(z) 



(40) 



From Eqs.(p7|) and ( J3q ) one notes that F cxc = 1 is the free energy cost of the interface in our units. 

Next we consider a semi-infinite system, in order to discuss wetting behavior in the strong segregation 
limit. Analogously to Eq.rtl]) , we choose a bare surface free energy of the form 



^smf - j(Pa - Pb) - h APA - Pb) 



(41) 
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where p A , p° B are the densities at the surface, and hi is a "surface field" in suitable units, A being the 
analog of the " extrapolation length" used in Eq. (|]) . 

We now note that the above profile, Eqs. (^) and ([lO]) is cutoff by the surface ^3|,|55|], and thus it 
is convenient to introduce the angle a = it — 2<po that corresponds to the values p A , p^reached in the 
surface plane. The bulk part of the excess free energy can be written as (note A = still holds) 

r°° / r\rh\ pir /2—a./2 I 

F c b x u c lk = jf dz L±j = J d0 sin20 = - (1 + cosa) (42) 

and combining Eqs. (^ and (fl2|) the total excess free energy becomes 

^cxc = \ + cos a Q - hi) + cos2 "A- (43) 
This free energy takes a minimum for 

cos a = ^ (hi- ■ (44) 
A second order wetting transition occurs for A > 0, at 

hl = \ + \ (45) 

while for 1/A < one has a first order wetting (as is actually observed in the simulation J57| of the 
model shown in Fig. 3). 

Finally we are now in the position to consider the case of thin films of finite thickness D , where now 
the constant of motion A is nonzero. The quantity that we wish to calculate is wo (D), defined from the 
inverse slope in the center of the profile 



w (D) = (2d PA /dz)- 1 = (2d^/dz)- 1 

z=0 



*=o Vl + A 



(46) 



where we have used that 



dp A dq\ dq A . M ,^d(j) 

— — = -r^ = 2q A —— = stn(2(j>)— (47) 
dz dz dz dz 

and we note that in the center of the profile <j> — tt/A since there pa — Pb by symmetry. From Eqs.(|34|) 
and (|36| ) we have deduced that 

(g) 2 = [A/4 + /»)] = [A + sin 2 (20)] /4, (48) 
and thus the counterpart of Eq. (p2|) for the profile becomes 



/4 



\ dz 

' -.Ei (ir/2 - a, 1 ) if A > 0, (49) 



v'TTA V VTTA 



= E X f 7r/2-arccos( ^==^1 ,VT+Aj , if A < 0, (50) 

E\ denoting the elliptic integral of the first kind. Denoting the angle 4>{z = 0) = </>o> and correspondingly 
a = tt — 20o, the excess free energy can then be written as 
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F, 



bulk 



--AD + 4 



--AD + 2 

4 



<t>0 



tt/4 



tt/4 



d0 
dl 



A + sin 2 (2^) 



4AD + yTTA£; 2 (vr/2 



for 



F 



bare 



-2h\ cos a + — cos 2 a 
A 



From these results one can show immediately for D 
that 



where A — ► and 



A w 16exp(-D) - 8aexp(-L>/2) 



(52) 

7r/2 hence a — * 0, 
(53) 



As anticipated above, A can be either positive ( if a < 2exp (—D/2) ) or negative ( if a > 2exp (—D/2) 
). In this limit the total free energy excess that must be minimized becomes 



F 



total 



const-4aexp(-D/2) + a 2 (/ii-2/A+l/2), a <2cxp(-£>/2) 



(54) 



while for a > 2exp(— D/2) an extra term 4Z?exp(— D/2) [a — 2exp(— D/2)] has to be added in Eq.(53). 
We are interested in the wetting case hi > h\. In that case, minimization of F*°* al with respect to a 
yields 



which yields in Eq. 



a = 2exp(-D/2)/(/n - 2/A+ 1/2), 
hi - (1/2 + 2/A) 



A = 16exp(-£>) 



hi + 1/2- 2/A 



(55) 



(56) 



Noting that h\ = 1/2 + 2/A and wq (oo) = 1 in our units, the final result for the reduction of the 
interfacial width can be cast in a form very similar to Eq . (pl|) , namely 



MP) 

wo(oo) 



1 + 16 



(hi - hl)exp(-D/w (oo)) 



hi 



-1/2 



(57) 



Also for the opposite limit -D — > 0, a— >7r/2an explicit analytical result is easily derived, since in this 
limit 



and 



D 

T 



D 



— 2hi cos a + — cos a 



D « 2w (| - aj 

( 7T 
2/11 (2 



2 1 



(f 



Minimization with respect to a yields 



7T 

— — CK ! 



hi 



hiD, 



D —f 



2 ~" 2/X + l/D 
and hence we obtain an equation analogous to the weak segregation limit, 

w (D) » hx/2. 



(58) 



(59) 



(60) 



(61) 



A > 



-\AD + E 2 (vr/2 - arccos (^=) , \/TTa) + %Ei (vr/2 - arccos (^=) , \/TTa) for A < 0, 

E2 being the elliptic integral of second kind. The surface part of the excess free energy, from 
Eq.([ll]),can be written as 



(51) 
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IV. NUMERICAL SELF-CONSISTENT FIELD CALCULATIONS FOR CONFINED 
INTERFACES OF POLYMER MIXTURES WITH FINITE CHAIN LENGTH 



We start the treatment by writing the partition function Z of a system of ha chains of type A and 
us chains of type B in the volume V with interactions £(pa,Pb) as a functional integral pS 



Z = 



n A \n B \ 



V[f a }V[f a ]V[f }V[f^]exp 



dV £(p A ,p B ) 



(62) 



where T>[r a ] stands for the functional integration of the coordinates of all the monomers of type A, T>[rp] 
the corresponding term for the monomers of type B, and "Pl^y, V\r^\ are the corresponding probability 
distributions of the chain conformations in the reference system. Since we compare the results of the 
SCF calculations with simulation data (i^] we use the bond fluctuation model in the athermal limit as 
a reference system. In the framework of this coarse grained lattice model each effective monomer blocks 
all 8 sites of an elementary cube of the lattice, and working at $ = 1/16 where half of the available 
sites are filled corresponds to a dense melt 21 1. In the following all lengths are measured in units of the 
lattice spacing. 

Defining density operators as 



^ HA N A 

few =&mi s ( r 



(63) 



a=l i= 1 



where the first sum runs over all A chains and the second sum over all monomers of the a'th chain, 
being at positions r a ^ , and similarly pb(t)- The normalized energy expression can be written as 



{PA - PB) 



(64) 



where the first term is identical to the second term of Eq.(j26]), while the second term represents a pairwise 
monomer-monomer interaction of finite range Iq (in practice we take Iq = 16/9, when we try to represent 
simulation results of the bond fluctuation model). Here q is an effective coordination number, and e a 
normalized energy between a pair of monomers (2qe — \, the Flory-Huggins interaction parameter) | pl| . 
While in the previous section we have considered the limit where Na = Nb = N — * oo and C — ¥ °° 
(incompressible melt of infinitely long chains), we here relax both these approximations simultaneously, 
allowing both N finite (in the numerical example shown in Fig. 3 we have chosen N = 32), and the 
compressibility C also is taken finite (</> = 1/16 corresponds to £ = 4.1, as discussed elsewhere pl|]). 

Eqs.(^||) - ( p4| ) is a very general formulation of the statistical mechanics of polymers, which we here 
drastically simplify in terms of a mean field approximation. The free energy density / then can be 
expressed as (T = — fc^TlnZ). 



*£fv = it ln pA + £ ln pB + v. 1 d3r 8 ' /M • ! ' 1J ] 



1 

v 



d 3 r{w A pA + wbPb} - tt" Ih^^a] - -£f~ ^-z B \w B } 

1\A l^B 



PB 



(65) 



with pa = 1 — pb 
chains, 



n ^y A , wa, wb being the effective fields, and za,zb the partition functions of single 



1 

V 



V[r a ]V[f a ]exp 



An 



a = A.B 



(66) 



The self-consistent fields follow from: 

wa — C (pa + Pb - 1) - qe 



1 1, 2 3 2 



(pa - Pb) ■ 



(67) 
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PA 



p A V T)z a 
N A z A Vw A 



1 Na 



I exp 



^M) n (rg,>) 



(68) 



where Af is a normalizing denominator, Af = ^2 a V w exp — ^(^a,*) j the sum over a is a sum 

over a representative sample of (Monte-Carlo-generated) polymer conformations (in practice 7 • 10 A 
chains and 7 • 10 6 B chains are used), and the probability V w of the confined polymer, in unnormalized 
form, is 





exp(±e w n v 



(69) 



Here, V w — applies if any monomer falls outside of the walls of the film, while otherwise the weight 
depends on the number n w of monomers experiencing the potential due to the wall e w : the + sign 
applies for A monomers at the left wall or B monomers at the right wall, while the - sign applies for B 
monomers at the left wall or A monomers at the right wall. 

In practice the unknown functions /3.4(f), /?s(f) that result from the solution of Eqs.(|67|) and ( |68| ) 
are found by choosing a Fourier decomposition into a set of basis functions appropriate for the chosen 
geometry, 



f k (z) = V2e 



irfcz 

IT 



fc = l,2,-- 



(70) 



The coefficients of this Fourier decomposition are found iteratively by the Newton-Raphson method. 
Note that the most difficult part of the present numerical SCF scheme actually is the summation over 
the sample of polymer conformations, which is done on a CRAY T3E multiprocessor machine. 

Fig. 4 shows typical profiles resulting from this method for N — 32 and three choices of e. While for 
e = 0.03 we have a strongly segregated case quite comparable to the corresponding Monte Carlo results, 
for e = 0.016 we have a profile that already develops the shape characteristic of the weak segregation 
case (Fig. 1). Fig. 5 shows the thickness dependence of w(D) resulting for this model at e - 0.03, 
comparing the SCF results with the corresponding Monte Carlo results from Werner et al flsfl . Note 
that these Monte Carlo results are not the full width of the profiles shown in Fig. 3 - which include also 
effects due to capillary wave broadening that cannot be taken into account by the above SCF treatment 
- but rather are constraint by recording the mean square interfacial width only on a lateral length scale 
B = 8 lattice spacings (the system is divided into a grid of B x B subblocks, and the local position of 
the interface in each subblock is recorded separately to obtain this local width). The proper choice of 
the value of this grid size B is discussed elsewhere Es,E9[. 



V. DISCUSSION AND CONCLUSIONS 



In this paper we have considered the change of the "intrinsic" width Wq(D) of an interface between 
coexisting phases confined in a thin film between "competing walls" a distance D apart. We have 
considered first two limiting cases which can both be treated by "square gradient" -type theories, namely 
the Cahn-Hilliard theory of a weakly segregated symmetric binary mixture (Sec. II) and the Helfand 
theory of a strongly segregated incompatible polymer mixture in the limit of infinite chain lengths 
(Sec. III). In both cases the reduction of the interfacial width w (D) due to confinement can be easily 
worked out, and a region where wo(D) varies approximately linearly with D occurs in both cases. 
Treating finite chain lengths for polymers, a numerical scheme has been used (Sec. IV) which also allows 
the treatment of cases intermediate between weak and strong segregation (Fig. 4). In this way, the 
different regimes proposed for the concentration profile in Fig. 1 quantitatively could be demonstrated 
explicitly. 

There is one important drawback of our treatment, however: while all variants of mean field theories 
(like those presented in the previous sections) readily yield "intrinsic" interfacial profiles, the latter 
are not well-defined in the framework of rigorous statistical mechanics, and consequently there is no 
unique way to define them either in a computer simulation nor in an experiment on real materials. 
The actual interfacial width w(D) observed in both simulations and experiments exhibits an additional 
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broadening due to fluctuations in the local position of the center h(x,y) of the interface away from 
its average position (h(x,y)) (= 0, in our choice of coordinate system, where the plane z — of the 
(x, y, z)-coordinate system is halfway in between the confining walls). 

Werner et al. [H have discussed the extent to which one can describe this broadening in terms of a 
convolution approximation, 

r+D/2 

Pa(z) = / dh^{z-h)P D (h), (71) 

J-D/2 

where p'f (z — h) is the intrinsic profile, for an interface centered at z = h, and Poih) describes the 
probability that a deviation h from the average value (h) = occurs for a film thickness D. It then was 
assumed that this probability distribution is a Gaussian, Po{h) = exp(— h 2 /2s 2 )/V2irs 2 , and hence one 
finds for the total mean square width 

w 2 (D) = W 2 (D) + ^ S 2 (D) (72) 

In order to calculate the additional broadening due to the interfacial position fluctuations, s 2 (D), an 
approximation was used where the interface is described by an effective interface Hamiltonian describing 
capillary waves in a harmonic potential: 



T-L cS {h) = J dxdy 



a{D) 



dh\ 2 /<% N " 



dx J \dy 



where a(D) is an effective interfacial stiffness (which converges to the interfacial stiffness a = cr(oo) of 
a free unconfined interfaces as D — > oo), a is a constant, and k _1 is a decay length which is of the order 
of the correlation length ^ in the weak segregation case, but of the order of wp(oo) = b/^6~x in the 
limiting case of strongly segregated polymer mixture with N — > oo. Using Eq. (|73|) , one then obtains 
that s 2 {D) oc D for large D, since 

2 1 C* IB d 1 k(D)D 

3 ( } ~ 27 J a(D)q 2 +aeM~ K D/2) ^ 8^{D) + °° nSt { ' 

Note that short wavelengths need to be cut off if they are less than some length B, in order that 
the integral in Eq.(^) converges. For self-consistency, all fluctuations on length scales smaller than B 
must be included in the intrinsic width wq(D). However, there is no obvious theoretical recipe that 
would uniquely define this cutoff B, and thus the separation of interfacial fluctuation into "intrinsic" 
and capillary wave" -type is somewhat arbitrary. In the simulations, B = 8 lattice spacings was chosen 
simply for the reason that then wo(D — > oo) agrees with the self-consistent filed calculations, and as 
Fig. 5 demonstrates, there is then fair agreement between the SCF prediction for w (D) and the Monte 
Carlo observation for all D. However, apart from this fact one has no reason for not chosing B = 7 
or B = 9, for instance. Also the analysis of the capillary wave spectrum in the confined geometry did 
reveal a pronounced dependence of cr(D) on D in the same regime where wq(D) differs appreciably from 
u>o(oo). As yet, an analytical approach to accurately predict the interface stiffening (cr(D) is enhanced 
for small D) due to confinement is lacking. It is conceivable that also the constants a and k of the 
effective interface potential Vn{h) — aexp(—nD/2)h 2 /2 are no true constants but also depend weakly 
on D. 

Particularly cumbersome is the theoretical understanding of the crossover between the weak segre- 
gation case and the strong segregation limit of a polymer mixture. Just as two lengths control the 
interfacial profile jl2], the length wq{oo) = b/y/6x shows up in the center of the profile, the radius of 
gyration R g = b^J N/6 in the wings, we expect two decay constants in the potential Vd(K), namely 

V D (h) = ^h 2 {aexp{-D/w (oo)) + a' cxp(-D/2R g )} (75) 

While the amplitude a 1 of the second term vanishes in the strong segregation limit, when rrib — > 1, it 
decays much slower with increasing D than the first term, for large N, and, hence the interplay between 
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these terms is subtle. We expect that a similar expression will interpolate between our weak segregation 
result for w (D) (Eq@), and the strong segregation result Eq.(|58[). Similarly,it is not clear whether the 
cutoff length B should be of the order of w;o(oo) or of the order of R g , in this case. Finally, we remind 
the reader that the exponential variation of Vo(h) with D is only appropriate for short range forces 
between the walls and the molecules of the mixture, not for the - physically more realistic - long range 
van der Waals forces. Thus, our treatment is a first step towards the resolution of a rather complex 
problem only. But there is clear experimental evidence that the effects discussed here are indeed 

practically relevant. Thus we hope that our study will stimulate further efforts to understanding this 
problem. 
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FIG. 1. 



Schematic order parameter profiles m(z) versus z, for a symmetrical phase separated mixture confined 
between two competing walls, such that the left wall prefers the i?-rich phase and the right wall prefers 
the j4-rich phase. The upper part shows a large thickness D of the thin film, the lower part a smaller 
thickness. The definition of the intrinsic thickness-dependent interfacial width wq(D) is indicated. For 
further explanations cf. text. 
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FIG. 2. 

a) Plot of the reduced interfacial widths wq(D)/wq(oo) versus -D/£b f° r -Woo = 2, 3, and 4, and two choices 
of A/£;,: A/£b = 1 (full curves) and A/£h = 2 (broken curves). 

b) Comparison of the exact numerical result for Wo(D)/wq(oo) versus D/^ for the choice M x = 2, 
A/£b = 1 with the approximation, Eq.(20) (broken curve). 
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FIG. 3. 

a) Order parameter profiles m(z) versus z for films of thickness D — 16, 32, 48, and 64 for the bond fluc- 
tuation model of a symmetrical polymer mixture (Na = Nb = N = 32, caa — (bb — —£ab = —ksTe 
with e = 0.03 corresponding to T/T c f, = 0.48,6^ = ±0.1 being a wall-monomer interaction of square well 
type and a range of 2 lattice spacings), using aixLxfl geometry with two L x L surfaces along which 
periodic boundary conditions act. All lengths are measured in units of the lattice spacing, b) Interfacial 
width w plotted versus D. Note that for D < 20 we have w oc D here, while the value u>o(oo) predicted 
by the self-consistent field theory (SCF) is shown as an arrow on the ordinate. From Werner el al p8[. 
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FIG. 4. 

Plot of the intrinsic order parameter profile m(z) vs. z, in units of the radius of gyration for chain 
length N = 32 and three choices of e. Curves result from the numerical self-consistent field scheme, 
Eqs.(||)-©. 
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FIG. 5. 

Intrinsic width w(D) plotted vs D/R g for polymers of chain length TV = 32, and energy parameters 
e = 0.03, e w = 0.1. The curve shows the result extracted from the numerical self-consistent field scheme { 
Eqs.(|63|)-(70)}, while squares are the Monte Carlo data of Werner et al. For extracting the intrinsic 
width from the simulations a lateral grid size B x B (with B = 8 lattice spacings) was used 
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